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2° ■ Abstract 

(N 

f^ , We present preliminary results on the computation of the QCD running cou- 

pling constant in the MOM scheme and Landau gauge with two flavours of 
^ dynamical Wilson quarks. Gluon momenta range up to about 7 GeV (/? = 

o ■ 5.6, 5.8 and 6.0) with a constant dynamical-quark mass. This range already 

-L . allows to exhibit some evidence for a sizable l//i^ correction to the asymptotic 

D , behaviour, as in the quenched approximation, although a fit without power 

'"^ ' corrections is still possible with a reasonable x^- Following the conclusions 

of our quenched study, we take into account l//i^ correction to the asymp- 
X ; totic behaviour. We find A^"^ = 264(27)MeV x [a-i(5.6, 0.1560) /2.19 GeV], 

C^ ' which leads to as{Mz) = 0.113(3)(4). The latter result has to be taken as a 

preliminary indication rather than a real prediction in view of the systematic 
errors still to be controlled. Still, being two sigmas below the experimental 
result makes it very encouraging. 
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The non-perturbative calculation of the running coupling constant of QCD is certainly 
a very important problem. In pure Yang-Mills it has been performed with several different 
methods, the most systematic ones using the Schrodinger functional |]I[, and the gluon Green 
functions |@-^. It is noticeable that the latter two methods, although very different, end up 
with perfectly compatible values for Aqcd- 

Of course the real challenge is to compute a^ with dynamical fermions. This task has 
been undertaken using NRQCD several years ago U[-M and, once extrapolated to Mz leads 



to rather satisfactory values for as{Mz)- Recently, the QCDSF-UKQCD collaboration ||T^ 
and the ALPHA one [^ have reported progress in determining a^ with two flavours using 
relativistic lattice QCD and nonperturbatively improved Wilson fermions. 

In this letter we will report our work consisting in applying the Green function method 
estimate [0-§] with non-improved Wilson dynamical quarks. The principle of the method is 
quite simple since it consists in following the steps which are standard in perturbative QCD 
in the momentum subtraction scheme. This gives immediately a nonperturbative estimate of 
the coupling constant at different scales. Its running can be confronted to the perturbative 
QCD expectation. We use the MOM renormalization scheme which corresponds to using 
an asymmetric subtraction point : p^ = p| = /i^,p2 = 0. This scheme proved to give rather 
good signals and, in spite of the zero momentum, no pathology has been seen. 

From our study of the pure Yang-Mills case 0-0 we have learnt two main lessons: one 
is that a study of the asymptotic behaviour of a^ needs a large energy window, since the 
value of Aqcd we are looking for depends on the weak logarithmic dependence of as on the 
energy scale fi, the second is that the l//i^ correction can be sizable up to a large energy. 

We aim at computing as, Aqcd and the power correction term with two flavours of 
dynamical quarks. This requires, as we shall see in more details, an exploration of the two- 
dimensional {go,msea) bare parameter space. To this goal we have run lattice simulations 
on several 16^ lattices. Notwithstanding the modest volumes of these lattices, we realised 
that some interesting physics can already be extracted. Furthermore, the exploration of the 
bare parameter space provides us with new data which might be of interest for unquenched 
studies by the lattice community. This legitimates in our opinion a progress report which is 
the aim of this paper. 

I. OUR STRATEGY 

We have computed in the Landau gauge the two-gluon and three-gluon Green functions 
leading to a nonperturbative calculation of a^'^**(/i) in the well defined MOM schemes [0,0. 
At energies /i above 2.6 GeV we will fit this function by 
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(/X') = «s,pert(/x') 1 + ^ , (1) 



where as,pert(/^^) is the perturbative running coupling constant computed to four loops from 
some fitted Aqcd, and as,pert c//i^ is a power correction which has proven, in the Nj = 
case, not to be negligible up to 10 GeV, and was eventually traced back to an OPE con- 
densate < A^A^^ >. The reason for choosing as in eq. (|l]) a non perturbative correction 
oc as,pert(/^^)//^^ instead of simply oc l//x^ is twofolds: 



(i) Theoretically, an OPE study |^ including a computation of the anomalous dimension of 
the coefficient oi < A^ > leads to an expected energy dependence close to as,pert(/^^)//^^; 

(ii) Practically in the quenched as well as unquenched case the fit with as,pert(Ai^)/At^ is 
much more stable for changes of the energy window than the fit with l//i^. The 
former stability will be illustrated in table 0. 

Interestingly, a^*^** is at the same time both the goal of our study and a very useful tool: 
from the lattice simulations one extracts the continuum as up to small lattice artifacts; 
lattice spacing ratios are then fitted to preserve the continuity of as{fi) for the whole set of 
data. 

This program is performed on hypercubic lattices in order to simplify the necessary ten- 
sorial analysis of the Green functions^. In the Nf = case we combined (3 = 6.0, 6.2, 6.4, 6.8 
quenched lattice simulations, i.e. a lattice spacing ranging from ~ 0.03 fm to 0.1 fm, in [^Q, 
allowing to reach momenta up to 10 GeV. 

With dynamical fermions the physics depends on two parameters, (3 which represents 
the bare coupling constant and Kgeg, representing the bare dynamical-quark mass. A wide 
energy window is reached by combining simulations with different lattice spacings and the 
same renormalised dynamical- quark mass expressed in physical units. The problem is of 
course that we do not know a priori for a given f3 which Kgea corresponds to one given 
renormalised dynamical-quark mass in physical units. This needs as mentioned above some 
exploration of the (/3, Kgea) parameter space to find one or several lines of equal dynamical 
masses. In view of the computational cost of such an exploration we have chosen to perform 
it on a small volume, 16^. 

We now would like to sketch our strategy to compute the lattice spacings and the renor- 
malised dynamical-quark masses in this exploratory stage on a 16^ volume. We proceed as 
follows. We start from a calibrating set of parameters (/3, /tsea) for which some published 
results yield the inverse lattice spacing a~^ computed from some hadronic quantity, for 
example the p meson mass. We then estimate a~^ for other values of (/3, Kgca) by matching^ 
the value of as{fi)- This uses as an assumption that we may neglect the dependence of as on 
the dynamical-quark mass, at least in the mass range under consideration. This assumption 
is not more arbitrary than any other calibration based on, for example, the physical p meson 
mass, which neglects the unknown dependence of the p meson mass on the dynamical-quark 
mass. 

Once we have estimated the lattice spacings for all our lattices with different sets (/3, Usea), 
we estimate rusea from the ratio df^A'^/P^ where A^ is the axial current and P^ the pseu- 



^ This does not allow to compute the p meson mass, which is better performed on lattices longer 
in time direction than in space. Some consequences will be discussed later. 

^Our matching procedure was proved to be succesful when applied to quenched data, where lattice 
spacings are well known P] . 



doscalar density, computed for a valence quarkQ with the same bare mass as the dynamical 

quark (K^al = ^sea)- 

At the end of this procedure we can fix with a reasonable accuracy a set of couples 
iP, i^sea) which includes our calibrating lattice and for which f3 the dynamical-quark mass 
remains constant in physical units when (3 is varied. This knowledge allows for a preliminary 
analysis of a^ with two flavours, of the resulting A^~ and power correction term, and finally 
of as{Mz)- This will be presented in this letter. 

Still we do not forget that finite volume effects may be large in such a small volume, 
that our dynamical-quark masses are large, etc. We will therefore briefly discuss sources 
of systematic uncertainties at the end of this letter. We are however now in a position to 
launch the calculations on a larger volume 24^* and/or with lighter masses and correct for 
the biases of our present results. 



II. SOME USEFUL PERTURBATIVE FORMULAE 

We now proceed to establish the conventions and to introduce the formulae that we will 
use in the following. as,pert(/^^) in eq. (|I|) stands for the perturbative running coupling 
constant expanded up to the fourth loop and verifying (in this section we write a instead of 
as,pert(/^^) to simplify the notations) 



d ^ / /^O 2 I /^i 3 I /^2 4 I Ps 



a = -l—a + —,a + -—a + — — « . (2) 



c/ln/i \2'K 47r2 647r3 1287r^ 

In all schemes 

2 IQ 

(3^ = ll--Nf [3,=hl--Nf, (3) 

while /32, /^s depend on the particular scheme. The values we need for the MOM scheme can 
be found in ref. |1^. The exact integration of eq. (0) to the third loop, with the standard 



boundary condition defining the A parameter [^, leads to [Q 



/3l 
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X exp < , ^- 



arctan, ^ '^ - arctan ^ 



(4) 



where A^^-* denotes the conventional two loops formula: 



^We call valence quarks the quarks which contribute to the current densities and propagate in the 
gauge field background. The latter depends on the sea quark mass. The theory is unitary only if 
valence and sea quarks have the same mass. However we will also make use of K^ai ¥" '^sea as an 
intermediate step. 



A(e)^ 
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(5) 



and A = 2(3q(32 — 4:(3f > in the MOM scheme which we use. If one only retains the first 
correction coming from the perturbative fourth loop, it can then be written 



A4100PS _ ^Sloops g^p 
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(6) 



In the previous formula, of course, the use of A, a and (3's stands for the A parame- 
ter, the running coupling constant and beta function coefficients in the particular MOM 
renormalisation scheme. From now on we will systematically convert A into Aj^ using [Q 



A 



MS 



A exp 



1 /70 
22 VY 



22 



9 



A^ 



(7) 



No analytical expression can exactly inverse neither three-loop eq. (Q) nor four-loop eq. 
(^). The following formula gives nevertheless an approximated solution to the inversion of 
the perturbative expansion of eq. (0): 



as,pert(/X ) = -5— 

Pot 
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87if3, log(t) 
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log(t) - 1 



(8) 



where t = log(/i^/A2). The exact numerical inversions of, for instance, eq. (^, can be easily 
obtained; but, of course, such an exact inversion and the approximated solution in eq. (^ 
should only differ by perturbative contributions of order higher than four loops. 

III. THE FIRST ITERATION 



A. Lattice spacings 

The lattice parameters which we have used for our simulations are displayed in table | 
together with our estimates of the pseudocritical k's, Kpc, defined in subsection |1I1C| , of the 
the lattice spacings and of the sea quark masses. 

We will not repeat the method used to extract a^ from Green functions as it is exactly 
similar to what was done in the pure Yang-Mills case. The MOM scheme uses the "asym- 
metric" three point Green function, i.e. with gluon squared momenta (0,/i^,/i^) where 



a 


'^soa 


Volume 


l^'pc 


a-i (GeV) 


m^ea (MeV) 


5.6 


0.1560 


16^ X 32 




2.19(8) 




5.6 


0.1575 


16^ X 32 




2.38(7) 




5.6 


0.1575 


243 X 40 


0.15927(5) 


2.51(6) 




5.6 


0.1580 


243 X 40 


0.15887(4) 


2.54(6) 




5.6 


0.1560 


24^ 


0.16053(3) 


2.19(8) 


164(7) 


5.6 


0.1560 


16^ 


0.16048(13) 


2.19(8) 


164(7) 


5.6 


0.1575 


16^ 


0.1593(1) 


2.42(9) 


79(3) 


5.8 


0.1500 


16^ 


0.15672(6) 


2.45(13) 


325(18) 


5.8 


0.1525 


16^ 


0.15555(12) 


2.76(7) 


173(4) 


5.8 


0.1535 


16^ 


0.15522(9) 


2.91(18) 


103(16) 


5.8 


0.1540 


16^ 


0.15499(6) 


3.13(13) 


64(5) 


6.0 


0.1480 


16^ 


0.15272(7) 


3.62(10) 


391(12) 


6.0 


0.1490 


16^ 


0.15262(7) 


3.73(13) 


308(12) 


6.0 


0.1500 


16^ 


0.15238(4) 


3.78(14) 


213(3) 


6.0 


0.1505 


16^ 


0.15240(5) 


3.84(15) 


169(8) 


6.0 


0.1510 


16^ 


0.15207(3) 


3.96(16) 


96(4) 



TABLE I. 



Data taken from literature and first iteration estimates from our runs. Upc is defined 



in subsection QIC. The dynamical-quark masses are renormalised in the MS scheme at 3 GeV. 



/^^ = n(27r/L)^, n being an integer 0. Let us nevertheless recall that the "asymmetric" 
three-point Green function turned out to be more convenient for our purpose than the sym- 
metric one (i.e. with gluon square momenta (/i^,/i^,/i^)): the high accuracy achieved in our 
quenched study 0-^ gave us some evidence that the "asymmetric signal" was less noisy 
than the symmetric one, while being as reliable (more momenta available and no observable 
infrared pathology due to the zero momentum). 



In table |, we give the full set of runs performed and the preliminary values obtained 
for a~^ and msea- In the first four rows, we also quote sets of values taken from literature 
1^,[15|. In particular, the value a~^ = 2.19(8) GeV for (5.6, 0.1560) is taken from |T^ and 
we will use it to calibrate all our runs. As already mentioned, ratios of lattice spacings 
result from imposing the continuity of ^^(/i) from different lattices, neglecting the expected 
small dependence of a^ on the dynamical mass rrisea.- Since the error on the calibrating 



^Any integer verifies n 
square theorem). 



E. 



=0,3 



n- 



for at least one set of integers ni,i = 0,3 (Lagrange's four- 



propagates trivially to Aqcd we will use 2.19 GeV without its error until eq. (|T5]); a 



discussion of these errors will follow eq. (15). Thus the errors quoted here for a ^ only stand 
for the ratios. 



B. Sea- quark masses 

Once the lattice spacings are estimated, we also need to compute arrisea.- To this aim we 
compute the propagators of valence quarks for several Kvai among which one with Kvai = ^sea 
in order to be able to deduce the mass of the sea quark from the estimated mass of the 
valence quark. 

This is done using the ratio 



^ 1 Eg ^5 (0)^0^0 (X,t) 

^ 2 E.-A(0)P5(x,t) ^ ' 

where P5 is the pseudoscalar density, and A^ the axial current. 

To estimate the ratio p in (j^) we have used two methods. The simplest consists in 
looking for a plateau of the ratio, the time derivative in the numerator being computed 
by a symmetrised discrete difference. However this "brute-force" estimate appeared to be 
affected by some strong 0{a) effects in several cases. 



The second method, which is a variant of the one proposed in |[T6|, fits on some time 
interval the < P5P5 > in the denominator by a cosh function and the < P^Aq > in the 
numerator by a sinh with the same "mass" term]]. The time derivative of the sinh in the 
numerator is then proportional to the cosh in the denominator and the ratio gives an estimate 
of the ratio p in (^. We have used the second method because it turned out to be more 
stable against the change of parameters (domain of the fit) and to provide a better continuity 
when K, is changed. 

The valence mass is given by 

aw^vai = ^P (10) 



For Za and Zp in the RI-MOM scheme, we have taken |T7[ Z^ = 0.77(1) and Zp = 0.54(1) 
i.e. Z^/(2Zp) ~ 0.71 at /z = 3 GeV (the value of Zp is derived from the Ward identity value 
of Zp/Zs HH)- The large Goldstone pole contribution stressed in ref. |19| is claimed]^ to be 



eliminated in this value of Zp. A more careful study of the renormalization constant will be 
performed soon. 



^It is not really a mass since the very short time interval considered does not allow to isolate the 
ground state. It nevertheless turns out that the data for 5 < t < 11 can be satisfactorily fitted 
respectively with a cosh and a sinh. 

^Indeed, in ref. [23| this method is applied and agrees with ref. |24|, obviously free of Goldstone 



boson contribution. Although the Goldstone boson question is not commented in |23] this seems 
to confirm the above-mentioned claim. 



From arrisea. and a~^ we extract the masses presented in table |. These masses are 
computed in the MS scheme (3 GeV); the conversion from RI-MOM to MS is obtained by 
using formulae involving the four-loop anomalous dimensions of the quark mass [^ . 

C. Pseudocritical k: Kpc 

For any parameter set (/3, Kgea), having computed the valence masses for several values 
of Kvai we extrapolate to a vanishing valence mass. We call "pseudocritical k", i^pdi^sea), the 
value of Kvai for which rTiyai = 0. The values of arrival as a function of l//€vai are perfectly 
compatible with linear fits. 

The pseudocritical k's as a function of l/^sea are also compatible with a linear fit except 
for one point at /3 = 6.0: Kgca = 0.1510. We did not succeed to understand the reason 
for this unusual behaviour and have for the moment withdrawn this point from our fit for 
j3 = 6.0. We call "critical k" {hc{P)) for one P the value of Kgea at which the extrapolated 

Kp(. IS equal to Kgga- f^pc\l^c) ^ ^c- 

Our results for Kc are the following: 

fi;c(5.6) = 0.158480(32), ^^5.8) = 0.154682(34), fi;c(6.0) = 0.152012(34) (11) 

At j3 = 5.6, we made the extrapolation using two runs performed by us, (^sea = 0.1575 
and Ksea. = 0.1560 (24^) ) and a third one at Kgea. = 0.1580 taken from [^. The value of 
Kc(5.6) we obtained is perfectly compatible with the one (kc(5.6) = 0.15846(5)) published by 
SESAM [0. Replacing our run on 24^ by the one on 16^ induces no significant difference. 



To our knowledge, the last two values are new. 

D. Some tests of finite volume eff"ects 

It is clear that a critical point in everything we report here is the risk that finite volume 
effect might spoil our results. In our mind the present work should be mainly a preparation 
for similar runs on larger volumes (24^) and we want to be sure that the information gathered 
on 16'^ is relevant enough to tune our parameters for a larger volume. We performed two 
checks with this purpose. 

The first one is the comparison of the two runs at /3 = 5.6, Ksea = 0.1560 reported in 
table |. It can be seen that there is no significant difference between the results for 16^ and 
24^. 

The second one relies on the idea that there could be some kind of first order phase tran- 
sition at very small volume, a deconfinement and/or chiral restoration transition. Chiral 
symmetry restoration has the effect of eliminating the Goldstone boson and thus of invali- 
dating the relation mp oc {nig + rriq) where mp is the lightest pseudoscalar meson "mass"^ 
and rriq {rriq) the (anti) quark mass. 

Our analysis has found empirically that all our lattice dataj^ can be fitted to a good 
accuracy by the following formula: 



^We included the run at /3 = 6.0, Kgea. = 0.1510 in this analysis to test in a different way its chiral 



rrii 



25msea + 



V 



(12) 



Here rUg = rrig is the dynamical-quark mass, and V the lattice volume, both expressed in 
physical units. We still take the quark mass renormalised in the MS scheme at 3 GeV. We 
then obtained 



B = 2.74(5) GeV, r = 1.41(5)GeV2 fm^ 
from a best fit with a x^/d.o.f. = 0.57 (see Fig. [1|). 



(13) 
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FIG. 1. Finite size effects on chiral behaviour. 

Eq. (|I^) shows a strong but smooth finite volume effectQ without any sign of a sudden 
change of regime, as would be the case with a first order phase transition. 

In the infinite volume limit we should recover the pseudoscalar mass mp^oo- Indeed we 
checked for (3 = 5.6 and Kgea = 0.1575 that mp^^ ^ 0.432(18) GeV^, in good agreement with 
SESAM [|T|]: m^ ~ 0.432(9) GeV^. Furthermore m|,„^ = 2.74(mg + rUg), which from the 
pion mass gives (m„ + md)/2 ^ 3.6 MeV and from the kaon mass uis — 90 MeV at 3 GeV. 
This compares fairly well to other lattice estimates. 



behaviour; however it did not exhibit any pathology here. 
^'^We cannot compare our finite volume correction in eq. (p!2D to existing theoretical studies of 



finite volume effects on ground state energies |21|: indeed, as already stressed, mp is not really a 
ground state energy. 



IV. SECOND ITERATION: FITTING Aqcd AND POWER CORRECTIONS 

A. Fitting Aqcd and 0(1/^^) coefficient 

Once we have an approximate estimate of the lattice spacings and dynamical masses, we 
now proceed with a combined fit of Ug on the line with approximatively constant dynamical- 
quark mass which goes through (3 = 5.6, Ksea = 0.1560: /3 = 5.8, Kgea = 0.1525 and P = 
6.0, Ksea = 0.1505. This allows to reach momenta as large as ~ 7.0 GeV, large enough to see 
the asymptotic behaviour, provided we take into account 0(1/ fi^) corrections. 

According to our ansatz (|l|), we need to fit simultaneously the lattice spacing ratios, 
and the parameters Aqcd and c (coefficient of as,pert//^^)- To fit the lattice spacings one 
needs some analytic function to interpolate between the measured points and to adjust its 
parameters simultaneously with the lattice spacings to the smallest x^- Two approaches are 
possible. 

One approach is to use the asymptotic four loops behaviour plus as,pert//U^ corrections 
as the analytic function. This would allow to reach both goals with one stroke but at the 
expense of eliminating about one half of the points at /5 = 5.6, Kgea = 0.1560 which are too 
low in energy to follow the asymptotic behaviour. 

The other approach proceeds in two steps as follows. To fit the lattice spacings we 
have used polynomials. At /5 = 5.6, Kgga = 0.1560 we have used both the 16^ and the 24^ 
lattices. A universal polynomial (Fig. ^ Matching of lattice spacings) fitting all points of 
the four lattice settings considered does indeed exist except for a few points which happen to 
correspond to n = (L/i)^/(47r^)<2 — 4 where L is the length of the lattice. We attribute this 
behaviour to a strong finite volume effect and exclude the points below some IR cutoff. 
Varying this IR cut from n > 2 to n > 4 leads to a variation in x^ from x^ /d.o.f. = 1.06 
to x^ /d.o.f. = 0.79. For lower IR cut-offs the x"^ increases dramatically, while for higher IR 
cut-offs too many points are excluded. The uncertainty induced by the choice of the cutoff 
is taken into account in the systematic error which affects the values quoted below. At this 
point, it should perhaps be emphasized again that this procedure was tested on quenched 
data, providing us with the generally admitted lattice spacing ratios. In this way we obtain: 



a-\5.8, 0.1525) = 2.85 ± .09 ± .04 x ° 'twcev^^^ ^eV 
0-1(6.0, 0.1505) = 3.92 ± .11 ± .07 x ""'g-'g-y'"^ ^eV 



2.19 GeV ^^' (-lA) 

0-^5.6,0.1560) r^^^T V / 



where the central value corresponds to a cut at n > 3, the first error is statistical and the 
second is systematic. From now on, we are going to use these ratios, correcting the first 
iteration estimates shown in table |I[ 

Once the lattice spacings have been estimated we perform a combined fit of Aj-j^~ and 
the coefficient c as defined in eq. ([^) with a^l^^t given by the r.h.s. of eq. (||). The result is 
plotted in Fig. ||, {Asymptotic fit of as). 

From the results in table || we conclude: 

Ar - 2«(27)q^M^MeV .^2.7(1.2) 
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1(5.6,0.1560) 
2.19 GeV 



GeV 



(15) 
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FIG. 2. Fits of a, obtained with four lattices. 



The same analysis using the formula in eq. (|l]) leads to 

252(10)MeV c = 1.0(l)GeV2 



A^=° 



MS 



(16) 



from our Nf = data [] §. 

If we use the first approach mentioned above, i.e. fitting from the beginning with the 
formula of eq. (|I]), the result are perfectly compatible with eqs. (JT^) and (]T5|). If we try the 
same procedure without power corrections, keeping an energy window ranging from 2.6 GeV, 
we can obtain a best fit with the following parameters (see FiG. |^): 



a-^(5.8, 0.1525) = 2.86(6) x 
0-1(6.0,0.1505) =4.08(9) x 



-1(5.6,0.1560) 

2.19 GeV 
-1(5.6,0.1560) 

2.19 GeV 



GeV 
GeV 



and 



MS 



N 



345(6) X^/d.o.f. = 0.96 



(17) 



(18) 



It is not surprising that A— ^ is larger when the fit does not include l//i^ corrections, since 
the parameters A^~ and c vary naturally a contrario: if one increases the other decreases. 



There is a possible contradiction between the acceptable x^ found in (|T8D and the fact that 
in eq. (|T5D the coefficient c is three standard deviations away from 0. This might be due to 
some correlations in the data; further study is needed to settle the origin of the discrepancy. 



i^If the fit for Nr = is performed according to a^'^**(/i^) 



Nf=0 



a. 



s,pert 



(/^^) + -% instead of eq. ([l|) 



the result is A^ 



~ 237(10)MeV 
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/imin (GeV) 


Ams (MeV) 


c (GeV2) 


xVd.o.f. 


2.6 


264(24) 


2.66(77) 


0.58 


3.1 


256(20) 


3.03(85) 


0.56 


3.6 


267(29) 


2.51(1.21) 


0.54 


4.1 


269(29) 


2.44(1.75) 


0.71 



TABLE II. Four-loop fit with power corrections, eq. (ffl), on varying energy windows (> /imin) 
The stabihty of the fit is fairly good. 



Let us nevertheless underline three facts that make us confident about the pertinence of 
incorporing power corrections: 

1. The OPE analysis still holds in the unquenched case, so we theoretically expect the 
presence of power corrections; we wonder about their magnitude. 

2. In the pure Yang-Mills case 0-§] where high accuracy computations were achieved, 
the purely perturbative formula could not provide us with a good fit, but the one 
including power corrections managed to. This "numerical evidence" is comforted by 
the good agreement of our value of A^~ and that of the ALPHA collaboration []5|,^. 

3. Provided we do control the systematic uncertainties, the current data already show a 
tendancy to discriminate in favour of the inclusion of non-perturbative terms. Indeed, 
in the case of a Htwithout power corrections, the lattice data are lying over the fitting 
curve for the lowest values of our energy window, and under the curve for the highest 
energy values. We expect this tendancy to become clearer with higher statistics. 

Of course, we are conscious that the present letter would benefit from a deeper discussion 
of systematic errors. In particular the use of non- improved dynamical quarks leads to 0{a) 
errors. However our lattice spacings are all rather small and these errors partly cancel in the 
ratios of lattice spacings. The dominant error is thus an overall 0{a) error on the calibrating 
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lattice spacing. This error propagates multiplicatively to the values of A^-^^" and c. We are 
not in a position at the moment to estimate in a reliable way this systematic error. A rough 



estimate can be obtained either by analogy with the quenched case, looking at Fig. 4 in ||2^ 
^1 



for a ~ 2.2 GeV or directly in the unquenched case from Fig. 2 in WM. A crude estimate 



Nf 



is 20%. This would give 50 MeV on Aj^ . The effect on as{Mz) will be a systematic error 
of ±.004. 

It is useful to notice that the 0(a) errors are O(amsea) or ©(oAqcd)- No errors O^a/i) 
are expected due to symmetry reasons: the hypercubic symmetry of the lattice [fiu — *■ —f^u 
for any z/) implies that the momentum dependent errors are 0(a^/i^), exactly as in the pure 
Yang-Mills case. We did not see any significant effect of the latter errors for the momenta 
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0.7 



0.6 



0.5 



s 



0.4 



0.3 



0.2 



• p=5.6 

■ p=5.8 

P=6.0 

— Perturbative fit 




FIG. 3. Purely perturbative fit of a^ 



that we have considered. In particular such errors should show up in fig. 
deviation from the global fit for the data with largest a^^ . 



as a systematic 



B. Estimating as{Mz) 

At an energy of the order of the Z meson mass the 0(l//i^) power correction becomes 
irrelevant. We will therefore only keep as pert? the perturbative part of a^ from our fit and 
extrapolate. We proceed as indicated in [l^. We start from an energy of 1.3 GeV, the MS 
charm mass which is taken as the charm threshold^. At such an energy we will extrapolate 
from our quenched and two-flavour results to three flavours. We then start evolving up with 
four flavours to the beauty threshold, 4.3 GeV, and then further up with 5 flavours to Mz- 

Applying eqs. (§-^ with the values of Aj^ in eqs. (|15D and ([T6|) and assuming 
0-^5.6, 0.1560) = 2.19 GeV we get in MS scheme 



a 



.Nf=0, 
s,pert> 



1.3) = 0.259(6), 



a. 



Nf=2 
s,pert 



;i.3) =0.306(20), 



a. 



Nf=3 
s,pert 



where the Nj = 0, 2 results come from direct lattice estimates in 
the Nj = 3 has been extrapolated from the two latter[^ 



(1.3) = 0.329(26) (19) 

P] and in this work, while 



^^We have preferred to follow the tradition here, although it is not clear to us why one should use 
the MS mass and not the pole mass, and why the threshold is at rfic and not 2mc where the charm 
loop dispersive contribution starts for the gluon propagator. 

^^We simply assume that the extrapolation to an odd number of flavors is legitimate, not knowing 
what to do better. 
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The evolution up to Mz (where non-perturbative corrections are neghgeable) and down 
to Mt- gives 

a,{Mz) = 0.113(3)(4) as,pert(M,) = 0.283(18)(37) (20) 

where the second error conies from the systematic error on the cahbrating lattice spacing. 

V. DISCUSSION AND CONCLUSION 

We should reemphasize that this is mainly a progress report. Most of the results reported 
here were performed on small volumes and with rather large quark masses. Our goal was to 
undertake a first exploration of the parameter space. It turned out that the results seem to 
make sense. The rather smooth junction of the a^ points from three different lattices show 
that overwhelming ultraviolet or infrared lattice artifacts are absent. 

The points from different lattices with identical momenta do coincide unless 
(Lyu)^/(47r^)<2. Suffering presumably from strong finite volume effects these points have 
been excluded from the global fits. The comparison at /5 = 5.6, Kgca = 0.1560 of the 16^ and 
the 24^ volumes are encouraging and should be extended to other sets (/5, /tsea)- The finite 
volume effect on masses seems to be well accounted for by eqs. (|12D , (|^) , and the good 
agreement of mp^oo with the estimate in f^, performed on a larger time interval, confirms 
this optimism. 

Our result for as{Mz) is about 2 standard deviations below the world average experi- 
mental as{Mz) = 0.119(2) ||TB|]. It is slightly larger, although compatible within errors, with 
the result|3 of ^■. a,(Mz) = 0.1076(20)(18). Older results using NRQCD were closer to 
experiment: af^Mz) = 0.1174(24) §, af)(Mz) = 0.118(17) §. Our resuh for a,,pert(M^) 
is also 2 cr's below the experimental value of 0.334(22) MeV |2^. However, the meaning 
of this comparison is unclear because we cannot take into account the non-perturbative 
contribution to a^^ at Mr. 

We consider the fact that our preliminary result is 2 a's below experiment as very encour- 
aging. We should stress that the error presented in eq. (pO]) corresponds to the statistical 
error and only to some systematic errors: mainly the choice of the fitting window and the 
calibration error. Other systematic effects should be systematically explored such as that of 
the dynamical-quark action and that of the mass of the dynamical-quark (ours are rather 
heavy). A calculation with a lighter dynamical-quark mass is in progress. As a final remark 
we would like to stress that our value for as{Mz) is strongly correlated to the rather large 
1/yU^ corrections that we find in our fit. Starting from eq. (|T^) i.e. from a fit without power 
corrections we obtain as{Mz) = 0.1211(3)(40). As already stated, in the fits, Aj^ and c 
show an understandable tendency to vary a contrario. We are clearly encouraged to follow 
on this analysis and try to refine our result for as{Mz). 



^'^This results from the fact that our value A^ = 264(27) is larger than the value 217(16) (11) 



from [|10|. 
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Three-Gluon Coupling with Parametrization alplia4_ope: 



0.65 



0.6 



0.55 



0.5 



0.45 



0.4 



0.35 



0.3 



0.25 



0.2 



A^''ms= 0.331(2.) 
P, = 0.00(0.) 
a^XO.O) = 3.94(0. 
Q^XO.S) = 2.80(0. 
a^YS.e) = 2.19(0. 
X" = 0.89 
Npfit = 40 




2.5 



3.5 



4.5 



5.5 



Semi-exact 4-loop MOM~ Scheme - ( 2.60 < |^ < 8.00) 



